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We solve the problem of the transmission and reflection of phonons and rotons at the interface 
between superfluid helium and a solid, for all angles of incidence and in both directions. A consistent 
solution of the problem is presented which allows us to rigorously describe the simultaneous creation 
of phonons, R~ , and R + rotons in helium by either a phonon from the solid or a helium quasiparticle 
incident on the interface. The interaction of all Hell quasiparticles with the interface, and their 
transmission, reflection and conversion into each other, is described in a unified way. The angles 
of propagation and the probabilities of creating quasiparticles are obtained for all cases. Andreev 
reflection of helium phonons and rotons is predicted. The energy flows through the interface due to 
phonons, R~ , and R + rotons are derived. The small contribution of the R~ rotons is due to the small 
probability of an R~ roton being created by a phonon in the solid, and vice versa. This explains 
the failure to directly create beams of R~ rotons prior to the experiments of Tucker and Wyatt in 
1999. New experiments for creating R~ rotons, by beams of high-energy phonons (h-phonons), are 
suggested. 

PACS numbers: 47.37.+q 



I. INTRODUCTION 

Many physical properties of continuous media at low 
temperatures can be described in terms of quasiparti- 
cles. The quasiparticles of superfluid helium are called 
phonons, R~ rotons, and R + rotons. They have a non- 
monotonic dispersion curve and the R~ rotons have a 
negative group velocity, i.e. their momentum is directed 
opposite to the group velocity, see Fig.l. The phonons 
and rotons are observed in many experiments, such as in 
neutron scattering in helium [1] and in the direct exper- 
iments [2, 3] where beams of superfluid helium quasipar- 
ticles are created by a heated solid. The quasiparticles 
propagate in the helium, and interact and reflect from 
different surfaces. Also they quantum evaporate helium 
atoms from the free surface. These have been investi- 
gated both experimentally and theoretically (see for ex- 
ample [4-8]). Interestingly, R~ rotons were not detected 
in direct experiments until 1999, when they were finally 
created by a specially constructed source [3]. They were 
observed by quantum evaporation. All the earlier at- 
tempts to create R~ , with ordinary solid heaters, were 
unsuccessful. 

The problem of the interaction of rotons in superfluid 
helium with interfaces, their reflection, transmission, and 
mode change, was first considered in [9]. However, the 
method used there did not take into account the simul- 
taneous creation of phonons, R + and R~ rotons by a 
phonon in the solid incident on the interface, and it 
could not distinguish between the R + and R~ rotons. 
Later in [10] it was shown that R~ rotons cannot be 
created at the interface with a solid by a phonon from 
the solid, provided we can neglect the possibility of cre- 
ation of the other quasiparticles in the same process at 
the interface. In the current work a consistent solution is 
introduced, which allows us to rigorously solve the prob- 
lem of the simultaneous creation of phonons and rotons. 
Also it describes in a unified way, the interaction of all 



Hell quasiparticles with the interface: their transmis- 
sion, reflection and conversion into each other. These are 
the fundamental elementary processes that determine the 
heat exchange between HeU and a solid, and the associ- 
ated phenomena, such as the Kapitza temperature jump 
(see for example [11]). We investigate all these phenom- 
ena. The probability of creation of each quasiparticle at 
the interface is derived for all cases. The failures of at- 
tempts to detect R~ rotons prior to experiments [3] is 
explained, and predictions are made for new experiments 
on the interaction of phonons and rotons with a solid and 
the creation of R~ rotons at the interface by high energy 
phonons (h-phonons) . 
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FIG. 1. The solid line is Q(k) from (1) for s = 230.7m/s, 
kg = 1.9828A" 1 , and A = -0.9667; the dashed line is the 
measured dispersion curve of superfluid helium [1] at the 
saturated vapour pressure. 



We describe superfluid helium with its distinctive dis- 
persion relation Q(k), with the maxon maximum and ro- 
ton minimum, within the framework of the theory de- 
veloped in [10]. The quantum fluid is considered as a 
continuous medium at all length scales. This model is 
based on the fact that the thermal de Broglie wavelength 
of a particle of a quantum fluid exceeds the average inter- 
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atomic separation. Then the variables of the continuous 
medium can only be assigned values, at each mathemat- 
ical point of space, in a probabilistic sense. 

The idea to describe superfluid helium as a continuous 
medium at microscopic scales has been successfully used 
for decades. Atkins [12] used it in the 1950s to describe 
the mobility of electrons and ions in Hell, when he intro- 
duced bubbles and snowballs of microscopic size. Lately 
the vortices in superfluid helium with cores of sizes of the 
scale of interactomic distances are being studied exten- 
sively, see for example Ref. [13] and the references cited 
there. Some recent simulations on the dynamics of atoms 
in helium nanodroplets [14, 15] also affirm that Hell is 
well described as a continuous medium at microscopic 
scales. 

As shown in [10], application of the methods of theory 
of continuous medium at microscopic scales demands the 
relations between the variables of continuous medium to 
become nonlocal. In the work [16], the nonlocal hydro- 
dynamics was introduced to describe small oscillations in 
superfluid helium, and in [17, 18] it was used to describe 
ripplon-roton hybridization and dispersion relation of rip- 
plons. The nonlocality allows one to analyse a continu- 
ous medium with an arbitrary dispersion relation. This 
posibility was discussed in [19]. However, the theoretical 
justification of this approach remained on the intuitive 
level until the work [10]. 

In this paper, following [10], the quasiparticles are de- 
scribed as wave packets propagating in the superfluid. 
The long wavelength excitations are phonons, while the 
short wavelength ones are rotons. R~ rotons correspond 
to the descending part of the dispersion curve and have 
negative group velocity, i.e. they propagate in the di- 
rection opposite to their momentum. This simple model 
allows us to use the methods of the theory of continuous 
medium and avoid the difficulties that appear in other 
phcnomenological models, such as [17, 18]. 

So, with the help of boundary conditions for the con- 
tinuous media at the interface, we find the creation prob- 
abilities of all the quasiparticles' creation when any of 
them is incident on the interface, as the energy reflec- 
tion and transmission coefficients for the corresponding 
wave packets. The method allows us to obtain the an- 
alytical expressions for the probabilities as functions of 
angles and frequency. 

In section 2 of the paper we formulate the problem 
and obtain the general solution of the nonlocal equations 
of the quantum fluid in the half-space. We consider a 
parameterised dispersion relation that is a good approxi- 
mation to the measured dispersion curve of superfluid he- 
lium. The solution is sought in the form that generalizes 
the solution for a monotonic dispersion of a general form 
obtained in [20]. The consequences of using the bound- 
ary conditions are discussed, these include multiple crit- 
ical angles, backward refraction and retro-reflection (or 
Andreev reflection [21]) of phonons and rotons (see also 
[18]). 

In section 3 the boundary conditions are used to de- 
rive both the amplitude and energy reflection and trans- 
mission coefficients for any incident wave for arbitrary 
incidence angles. The preliminary results for rotons at 
normal incidence were discussed in the authors' report 
at the conference [22]. 

In section 4 the energy flows through the interface due 
to phonons, R~~ , and R + rotons are calculated as func- 
tions of temperature. The contribution of the R~ rotons 



to the energy flows, in both directions, are shown to be 
small. This means that R~~ rotons are hardly created by 
a solid heater and are poorly detected by a solid bolome- 
ter. This explains why R~ rotons could not be detected 
in direct experiments until the work [3] . There they were 
created by a source made up of two heaters facing each 
other which allowed mode changes, and detection was 
achieved by quantum evaporation. 

The results obtained in this work can be used in other 
fields of physics. In particular they are important for 
classical acoustics, where the problem of wave transmis- 
sion through an interface has been solved only for the 
case when the dispersion relations of both adjacent media 
are strictly linear (see, for example, [23]), let alone non- 
monotonic. This problem concerning real media, with 
nonlinear dispersion, was of interest in the middle of the 
last century [24] and is still relevant today [25, 26]. 



II. DERIVATION OF EQUATIONS AND THEIR 
SOLUTION 

A. Problem Formulation 

Let us consider two continuous media separated by a 
sharp interface z — 0. In the region z < there is an 
ordinary continuous medium with sound velocity s so i and 
equilibrium density p so i. For the solid we only take into 
account longitudinal waves. 

Taking into account the transverse waves can be done 
in the same framework of theory of continuous medium 
and does not present any difficulty. However, the cal- 
culations become much more cumbersome, while on the 
whole the situation does not change. Due to the very 
small impedance of the solid-helium interface (see section 
2.4 and below), the reflection coefficients hardly change 
at all. For the transmitted waves additional critical an- 
gles appear, corresponding to the sound velocity of the 
transverse waves. Also it should be noted that taking 
into account both the longitudinal and transverse waves 
in the solid allows one to consider the contribution of 
Rayleigh waves, which give contributions to the trans- 
mission coefficients of He II quasiparticles into the solid 
at fixed incidence angles. For phonons with linear dis- 
persion this problem was solved in [27].. This problem 
for the helium-solid interface may be the subject of next 
paper. 

The region z > is filled with the quantum fluid with 
equilibrium density po and dispersion relation 51(fc) such 
that 

^) = ^ 2 {i + 2a| + !}. (1) 

Here s is sound velocity at zero frequency, k g is wave 
vector that determines the scale of the curve and A de- 
termines the form of the curve. For a range of param- 
eters, this relation is a good approximation to the mea- 
sured non-monotonic dispersion relation of superfluid he- 
lium (see Fig.l). For A < — 1 there are real k such that 
fJ 2 (fc) < 0, which is nonphysical, and for A > — a/3/2 
the curve is monotonic. For A 6 (—1,-^/2) the curve 
has the roton minimum at k = k rot and maxon maxi- 
mum at k = kfyidx 

as it should. We adopt the follow- 
ing set of values: s = 230.7 m/s, k g — 1.9828 A , and 
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A = —0.9667. Then the dispersion curve has the follow- 
ing parameters: the coordinates of roton minimum k ro t = 

0.9670kg = 1.913 A" 1 and ^ rot = Ml(k rot )/k B = 8.712 K 
is the Boltzmann constant); the maxon maximum is 
^max — f^(kmax) / ks = 13.8 K. These values are the ex- 
perimentally measured parameters of supcrfluid helium 
dispersion curve at the saturated vapor pressure [1]. 

We describe this quantum fluid by nonlocal hydrody- 
namics as developed in Ref. [10]. Accordingly, the quan- 
tum fluid, as well as the ordinary fluid on the other side 
of the interface, obeys the linearized equations of contin- 
uous media 



dp 

at 



-poVv 



~dt 



Pa 



-VP, 



(2) 



where v is hydrodynamic velocity, p and P are the devi- 
ations of density and pressure from the respective equi- 
librium values (for brevity we refer to them below as just 
density and pressure). The difference is that pressure 
and density in the quantum fluid are related through the 
non-local relation 



p(r)= J dV/i(|r-r'|)P(r'), 



(3) 



z'>0 



in which the integration domain is the region filled by the 
quantum fluid [28] . 

The suggested model well describes the interface be- 
tween superfluid helium and a solid, because for solids 
the relationships are local. In the frequency range of 
the dispersion curve of Hell, the dispersion laws of most 
solids, such as the heater materials of copper or gold, are 
very close to linear and they can be described as ordinary 
continuous media. 

Equations (3) and (2) lead to the integro-differential 
equation for pressure 



AP(r,t)= J d s r'h(\r-r'\)P(r',t), 

z'>0 



(4) 



that is set for x,y,t € (—00,00), z e (0,+oo). In the 
infinite medium, when the integration and definition do- 
mains are infinite, the Fourier transform of (4) gives us 
the relation between the Fourier transform of the kernel 
h(r) and the dispersion relation of the fluid f2(fc) [10]: 



h{k) = 



n 2 {k)' 



(5) 



For the dispersion relation (1) we obtain from the Fourier 
transform of Eq. (5): 



h(r) 



k, 



1 



47rs 2 r fci — k 2 _ 



( e ik + r — e~ ik ~ r ) 



(0) 



where k + and (— fc_) are the poles of h{k) in the upper 
half-plane C + 

fc± = k g (Vl - A ± iVl + A) /a/2, k + = k*_ <= C+. 

(7) 

Here the asterisk denotes the complex conjugate and C + 
is the upper complex half-plane. Due to the last condi- 
tion, the kernel, Eq. (6), despite the complex notation, 
is real. 



There is no convolution product in Eq. (4), either 
in the sense of one- or two-sided Fourier transform or 
Laplace transform, because the lower limit by z' is finite 
while the kernel is symmetrical h(r) = h(r). 

We consider the problem of waves transfering through 
the interface. As the equations (2) are local and coincide 
with the notation used in the equations of ordinary ideal 
continuous medium, the two boundary conditions on the 
interface (local!) are obtained from their integral forms, 
in the usual way, using the theory of continuous medium: 



P(x,y,z = -0,t) = P(x,y,z = +0,t), 
V x (x,y,z = -0,t) = V z {x,y,z = +0,t). 



(8) 



By applying the solutions of the equations of continuous 
media, on the both sides of the interface, the boundary 
conditions (8) give us the solution in the whole space and 
thus provide us with all the coefficients of reflection and 
transmission. The solution in the solid is well-known, 
and the solution in the quantum fluid is derived in the 
following subsection. 



B. Solution of Eq. (4) in half-space 



The equation that determines the relation between k 
and oj 



n 2 {k) 



(9) 



with Q 2 (k) from Eq. (1) is sixth order with respect to 
k. Its six roots are functions of oj and are denoted k^ 
for p = 1, . . . , 6. We note that if we used more terms in 
the polynomial f2 2 (fc), then the higher order equation for 
k would give 6 real roots and the other roots would be 
imaginary. 

In the problem of waves transfering through the inter- 
face, the two boundary conditions (8) can be satisfied for 
all x, y, and t only if all of the waves present on both sides 
of the interface have the same frequency oj and tranverse 
component of wave vector k T . A single monochromatic 
wave is not a solution of Eq. (4). Therefore, as there are 
in total six roots of Eq. (9), we search for the solution 
as a sum of six monochromatic waves, with the same fre- 
quency oj and transverse component of wave vector k T 
(the y axis is chosen along k T ), i.e. with the form 



P(r, t) = ^2A t _ l exp [«(k M r - 1 



,t)} 



Here the vectors are 



k 2 - k 2 



kl 



(10) 



(11) 

(12) 



The transverse component k T is real for physical reasons, 
but k^ and k^ z can be either real or complex, as we have 
to ensure the boundcdness of our solution only in the 
half-space z>0. 

After substitution of Eq. (10) into Eq. (4), we obtain 
the system of equations for the amplitudes and the 
equations for fc„ z , as functions of oj and k T . The system 
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for Afj, is 



E MK* - Kz)- 1 = o, 
E ^(*,« + fc-»)- 1 = o ) 



(13) 



where 



h 2 - k 2 

K ±z — K ± 



k T , fc+z — k_ z € C H 



The equations for k^ z are reduced to the form 
& (k 2 = k 2 + k 2 ) = lu 2 for /* = !,..., 6. 



(14) 



(15) 



The system of two homogeneous equations for the am- 
plitudes Afj, (13) ensure that no (nontrivial) solutions ex- 
ist with less than 3 non-zero amplitudes A^, i.e. there 
are no eigensolutions of Eq. (4), in the half-space, con- 
sisting of less than 3 monochromatic waves. This is the 
consequence of the nonlocality, which changes equation 
(4) itself in the presence of the interface. In infinite space, 
on the contrary, the domain of integration is the whole 
space, and the equation is solved by the Fourier trans- 
form, and its solution is a superposition of plane waves 
with dispersion (1). 

A solution of (4) , with the smallest possible number of 
waves being three, is constructed in the form (10), with 
three terms out of six, by picking a subset of any three dif- 
ferent roots {k a ,kj3, fc 7 } out of the set of six {fc M } /i= i i . ..,6- 
Then it can be rewritten with the help of Eqs. (13) in 
the form that contains a single amplitude: 

>(Q) 



P{k a ,k ff: k-,}(r,t) - P a/ 3 7 ~ 

[kaz kj rZ )(k az + fc_ 3 
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k/3z){k a z 


k~f Z ) 


{kp z 


- k +z )(kp z 


+ *-,) e , 


(kp z 


~ ky Z )(k/3 Z 


— k az ) 


(k^z 


~ k+ z ){kyz 


+ k. z ) r 



ik az z 



lk ^ z + (16) 



where Pf^ is chosen so that P{ ka ,k^k 1 }(r = 0,t = 0) 

p(0) 
r af3-f- 

The velocity is obtained from Eqs. (16) and (2): 



V{fe„,fc^,fe T }(r,i) 
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a/37 
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+ k- z ) 


UJ 


(k^z 


kaz ) (k^fz 


- kf3z) 



e ik{, z z_^_ 



(17) 



J(k T y-uit) 



As the two conditions (13) restrict the number of free 
amplitudes in Eq. (10) from six to four, any four linear- 
independent solutions of the form (16) constitute the ba- 
sis set of solutions of Eq. (4) for given u and k T , and any 
solution consisting of 4, 5, or 6 monochromatic waves can 
be represented as their linear combination. 



Roots of dispersion equation, 
out-solutions 



In- and 



The roots of Eq. (9) with il 2 (k) from Eq. (1) with 
respect to k 2 are k 2 = k 2 ^i for i = l,2,3, where are the 
three dimensionless roots of the cubic equation 



e 3 + 2A£ 2 +£- X 2 



0. 



(18) 



Here x = w/(sfc ff ) is the dimensionless frequency; £j(A, x) 
are some elaborate complex-valued functions. 

The most interesting frequency range is % € 
(Xrot,Xmax), where Xrot(A) and 

X.max 

(A) are the dimen- 
sionless frequencies that correspond to the roton mini- 
mum and maxon maximum respectively. For such fre- 
quencies there are three types of running waves in the 
quantum fluid, corresponding to phonons, i£~, and R + 
rotons. The branches are numbered in this case in the 
ascending order of the absolute values of their wave vec- 
tors kf. < k\ < kmax < ki < k rot < ks, so that i = 1 
corresponds to phonons, i = 2 to Rr rotons, and i = 3 to 
R + rotons. 

We now consider the problem of quasiparticles transfer 
through the interface. The quasiparticles are treated as 
wave packets that propagate in the two media. There- 
fore, when we build the solutions in the quantum fluid, 
we have to take into account that wave packets, as well 
as quasiparticles, propagate with their group velocities 
dfl/dk [29]. So, a wave packet of the quantum fluid com- 
posed of waves with wave vectors close to ko, with its 
length fco < kmax (so that it is a phonon wave packet) 
and the zth component kg z > 0, propagates away from 
the interface; but a wave packet, composed of waves with 
wave vectors close to one with length ko € (k max , k rot ) (so 
it is an R~ roton packet) and the zth component k$ z > 0, 
propagates towards the interface. 



Basis set of four solutions: 

phR R* phR R* phR"R* 



ph R~ R + 



P-lii. P^-iil, PL 2) : Ul . p «n : til 

Phonon incident from the solid: 

Po.= til 
Excitation incident from the helium: 



phonon: 


P t+ P (1> = 

1 out 1 in 


ph. 


R~ roton: 


P t +P (2, = 

1 out 1 in 




R + roton: 


P t +P (3, = 

1 out 1 in 





FIG. 2. The basis set of solutions in the superfluid 
helium in the half-space and their sums, which correspond 
to the different incident excitations. 



Let us construct the solution in the quantum fluid P ou t 
(the "out-solution") that is realized when a wave in the 
solid is incident on the interface. This solution should 
contain only such waves that constitute wave packets 
traveling away from the interface (i.e. waves with posi- 
tive group velocity) or waves that are damped at z->+oo. 
Picking 3 out of 6 vectors k^ is the same as picking their 
normal components k^ z , as they all have the same k T . 
The six normal components k„ z , obtained as solutions of 
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Eq. (15), are grouped into three pairs of roots ± y/kf — fc 2 . 
for i = 1,2,3. For each pair with the same i, either 
both roots are real which occurs for small enough an- 
gles k T < hi, or both roots are imaginary for k T > fcj. 
In the first case one root corresponds to a wave trav- 
eling towards the interface, the other to a wave trav- 
eling away from it. The second case, one root gives a 
damped wave in z > 0, while the other gives an exponen- 
tially unbounded wave. So, P ou t contains no more than 
three waves (and no less because there are no such solu- 
tions) and therefore has the form of Eq. (16) (see Fig. 2). 
The squared normal components of the three constituent 
waves are 

kf z = kf~k 2 . fori = l,2,3. (19) 

We define the signs of roots k iz for P out to be made up 
of waves with the normal components of wave vectors 
equal to k\ z , k 2z and k 3z . Then taking into account the 
negative group velocity of R~ rotons [29] , for the signs of 
real k iz we obtain fc lz , fc 3z > and k 2z < 0. If k T > ki for 
some i, the corresponding wave ~ exp (iki Z z) is bounded 
in z > for Im ki z > 0. Then in the general case we have 

fc r e(0,fci): 0<k lz <(~k 2z )<k 3z 
fcr : 0<{-k 2z )<k 3zi k lz eC+ (20) 
k r €(k 2 ,k 3 ): 0<fc 32 , k lz ,k 2z eC + . 

Then the out-solution has the form of Eq. (16) with the 
three wave vectors picked from the set of six with normal 
components k\ z , k 2z and k 3z : 

Pout = P{k a ,kf3,k^} [k az = ki Z7 kp z = k 2z , k lz = k 3z ] . (21) 

In order to solve the problem of a wave transfering 
through the interface from superfluid helium into the 
solid, we need also solutions containing waves that are 
traveling towards the interface (i.e. wave packets com- 
prised of these waves should be traveling towards the in- 
terface). We define them in the way, which is illustrated 

by Fig. 2. The solution P^J is constructed of waves with 
zth components of wave vectors (— ki z ), k 2z and k 3zi with 

(2) 

the amplitudes related through Eqs. (13). Solution P> n ' 
is constructed of waves with ki z , (—k 2z ) and k 3z . The 

(3) 

last one, P in contains waves with k\ z , k 2z and (—fez). 
Then the three sorts of in-solutions, that correspond to 
the three types of the incident waves, can be written in 
the form 

P$ = Pout\ kuM - kig) fari = l,2,3. (22) 
The Pj^ solution corresponds to the incident wave of 

(2) 

type i. So, a linear combination of P ut and P> n consists 
of one R~ roton wave (i = 2) that corresponds to the R~ 
roton wave packet incident on the interface, and all three 
waves that correspond to the reflected phonon, R~ , and 
R + roton wave packets. The amplitudes of the phonon 
and R + roton waves are the sums of the amplitudes of 

(2) 

those waves present in both P out and P> n ' . This is the 
solution in z > realized when an R~ roton is incident 
on the interface. The four solutions (21) and (22) are 
linearly-independent due to their structure and can be 
used as the basis set of solutions as mentioned at the end 
of the previous subsection. 

For x£ (Xrot,Xmax) the roots k lz are fully defined by 



Eqs. (19) and (20). For x& (0, Xrot) the roots £2,3 are 
complex and £2 = £3; then k 2z and k 3z are defined so 
that the roton waves ~ exp (ik 2z ^ 3z z) are damped, so 
k 2z = —k$ z G C+. In the limit x ~ ► the phonon waves 
have almost linear dispersion k\ w w/s, and it can be 
shown that k 2z — > k +z and k 3z — > — fc_ 2 . Therefore the 
amplitudes of all roton waves, that contain multipliers 
(k 2z — k +z ) and (k 3z + k- z ), go to zero, see (16), and 
the general solution tends to the ordinary superposition 
of incident and reflected phonon waves, with wavelengths 
much greater than the scale of nonlocality Ifcil" 1 . This 
limiting case can also be obtained by passing to the long- 
wave limit h(r) — > 8(r)/s 2 in Eq. (4), thus making it a 
local wave equation. 

The equation (4) was solved earlier in work [20] for the 
case of arbitrary but monotonic dispersion relation. The 
Wiener and Hopf method was used there, the application 
to this type of problem was suggested in [28] and devel- 
oped in [30]. It is much more general and seems to be 
more rigorous than the method used here, though less 
straight-forward. The solution of [20] can be generalized 
to the nonmonotonic case and can be shown to yield for 
the dispersion relation (1) exactly the solutions (21) and 
(22). 

In the work [17] the problem analogous to (4) was con- 
sidered in order to investigate the hybridization of ro- 
tons and ripplons. There the problem in half-space was 
replaced by the one in the infinite medium with sym- 
metrical extrapolation of the solutions with respect to 
variable z. However, in contrast to the usual differen- 
tial (local) equations, for an integral equation such as (4) 
the solution cannot be thus extrapolated symmetrically 
to zjO. Indeed, if we consider formally the solution of Eq. 
(4) in z < 0, it would be defined unambiguously by the 
solution in z>0 (from Eq. (16)), through the integral of 
Eq. (4). Direct substitution shows that the full solution 
is not even. This might be the reason that the method 
[17] gave wrong results near the surface waves threshold 
and was then rejected; in the next paper by the authors 
[18] another approach was used for that problem. 

D. Multiple critical angles and Andreev reflection 

Even before applying the boundary conditions (8) and 
finding the solution in the whole space, we can use the 
fact that two linear boundary conditions for the variables 
of continuous media are satisfied on the interface and 
derive a number of important consequences. First of all, 
we see that the solution is always constructed in a such 
a way that there are in total four outgoing (i.e. reflected 
and transmitted) waves; one in the solid and three in the 
quantum fluid. This is because the four conditions on the 
waves amplitudes, two from Eq. (8) and two from Eq. 
(13), can be all satisfied only when there are at least four 
outgoing waves. On the other hand, they can be at most 
four because, either in the formulation of the problem 
there is only one incident wave, or the requirement that 
the solution is bounded when some of ki Z are complex. 

Furthermore, the two boundary conditions (8) imply 
that all of the waves constituting the full solution have 
the same frequency u> and tangential component of wave 
vector k T . When wave i is propagating at angle 9i to the 
normal to the interface, k T =fcj sin 9i. Then if one of the 
waves is incident, the corresponding angle and k T are set 
and all the other transmission and reflection angles are 
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FIG. 3. When a phonon in the solid is incident on the 
interface, three quasiparticles, a phonon, R~ roton, and 
R + roton, are created in superfluid helium with the same 
uj and k T . The created R~ roton propagates backward in 
the transverse direction (i.e. retro-refracted). 



determined by the generalization of Snell's law 
sin 9 so i _ sinfli _ sin 9 2 _ sin 9 3 

Ssoi si s 2 S3 

Here s$ = oj/ki(uj) are the phase velocities of the cor- 
responding waves, that depend on frequency; k so / = 
ksoiz^z + k T e y is the wave vector of the wave in the solid 
and s ad = oj / k so i(uj) = const. The reflection angle for the 
wave of the same type as the incident one is equal to the 
incidence angle. 

From now on we will consider s so i > s, as is the case 
when superfluid helium is adjacent to a solid. Usually 
even the strong inequality holds. Then we have 

Ssoi > si > s 2 > s 3 > 0. (24) 

If a wave from the solid is incident at 9 so i, it is reflected 
at the same angle and the three waves are transferred 
into helium at angles 9i < 9 so i- The R~~ roton wave, as 
opposed to the others, due to its negative group velocity, 
propagates backward in the tangential direction (i.e. in 
the direction y— >— oo, see Fig. 3). 

Assume a wave i is incident from helium at 9i. Then 
according to Eq. (23) the transmitted wave in the solid 
has sin 9 so i < 1 for sin 9i < Si/s so i, and if the incidence 
angle is greater than the critical value, k so i z is imaginary 
and the wave in the solid is exponentially damped. Thus 
we obtain the three angles of full internal reflection 

sm6? = Si/s so i for « = 1, 2, 3. (25) 

In the same way there are three new critical angles de- 
fined for i>j (so that Si<sj): 

sin 6% = s t / Sj < 1 for {i,j} = {2, 1}, {3, 2}, {3, 1}. 

(26) 

If a wave i is incident and 9i < 0?J, then wave j (with j < i) 
has 6j G (0t,7r/2) and k jz G R. For 9, > 9% k 2 JZ < 0, the 
j-th wave is damped and the corresponding quasiparticle 
is not created. 

In an ordinary fluid, the group velocity of a wave packet 
is the same as the sound velocity and is constant. When 
a wave is incident on the interface the reflected wave has 



the same wave number and due to preservation of the 
transverse component of the wave vector k T , it is re- 
flected forward. This qualitative picture is maintained 
in the majority of all known physical systems. However, 
when the reflected wave (or quasiparticle) is qualitatively 
different from the incident, another possibility can be re- 
alized. So, when an electron of a normal metal is incident 
on the interface with a superconductor, it can be retro- 
reflected and converted into the hole of negative effective 
mass that travels back along the same line as the incident 
electron. This effect was discovered by Andreev (see ref. 
[21]) and is called Andreev reflection or retro-reflection. 

In our case when a helium quasiparticle is incident on 
the interface, three different quasiparticles are created in 
helium, with corresponding probabilities. If for example 
an R~ roton is created on the interface, along with other 
quasiparticles, when a phonon is incident, for the sake of 
brevity we will refer to this as "the phonon is reflected 
into R~ roton". While phonons and R + rotons behave 
in these processes like ordinary quasiparticles, R~~ rotons 
propagate in the direction opposite to their wave vector, 
due to the negative group velocity as already mentioned 
above. Therefore when a phonon or R + roton is incident 
on the interface, the phonons and R + rotons are reflected 
forward, while R~ rotons are reflected backward, or retro- 
reflected. In this way the transverse components of wave 
vectors k T are all equal. Likewise, when an R~ roton is 
incident, the phonon and R + roton are retro-reflected. 
Thus we have described the effect of Andreev reflection 
of helium phonons and rotons. 

If a monochromatic beam of phonons and rotons is 
incident on the interface at some angle 6, there will be 
up to seven reflected beams (see Fig. 4). We denote them 
as ij, which means "the beam of quasiparticles of type j 
created on the interface by the incident quasiparticles of 
type i" . The beams 11, 22 and 33 are reflected forward 
at the incidence angle 9, and therefore constitute a single 
beam ii. Beams 13 and 31 are also reflected forward, 
while 12, 21, 32 and 23 are reflected backward. Due to 
the relations (23), the beams 32 and 21 are reflected at 
angles greater than 9, and beams 12 and 23 at angles less 
than 9. 

The reflection angles depend on Si, which are functions 
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FIG. 4. When a single beam of quasiparticles is inci- 
dent on the interface, a set of reflected beams is created. 
The beam marked ij consists of quasiparticles of type j 
created by incident quasiparticles of type i. The beams 
ii are reflected specularly, while all others propagate in 
different directions. The beams 12, 21, 23, and 32 are 
retro-reflected. 
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of frequency. Therefore with an incident beam which 
is non-monochromatic, the reflected beams all become 
angularly diffused [20], except for the beam ii. However, 
the relations (24) hold at all frequencies, and therefore 
the qualitative picture is not modified. If we place a 
detector on the same side from the normal as the source 
at greater angles, it should register the R~ rotons of beam 
32 and phonons of beam 21, which were retro-reflected. 
Likewise the detector at smaller angles should register the 
R~ rotons of beam 12 and R + rotons of beam 23. Such 
an experiment could be carried out in order to verify 
qualitatively the current theory. 

A successful experiment would very much depend on 
the intensities of the beams to be detected, and therefore 
on the different creation probabilities for the quasiparti- 
cles at the interface. The derivation of these probabilities 
is the subject of the next section. 



III. REFLECTION AND TRANSMISSION 
COEFFICIENTS 



A. Phonon in the solid incident on the interface 

When a phonon in the solid is incident on the interface, 
it is reflected and three quasiparticles of different types 
are created in helium which travel away from the interface 
(see Fig. 3). The probability of quasiparticle creation is 
the fraction of incident energy which is reflected or trans- 
mitted as the corresponding wave packet. Amplitude re- 
flection and transmission coefficients can be derived in 
the approximation of plane waves [30]. 

In this approximation we consider a plane wave with 
frequency to and wave number k so i(ui) =U)/s so i incident 
on the interface at angle 9 so i to the normal. Then the 
solution in the solid is the sum of the incident and re- 
flected waves, the solution in the quantum fluid is P out 
from Eq. (21), which consists of three waves. All the 
waves have the same frequency w and transverse com- 
ponent of wave vector k T = k so i cos SO /. The pressure 
amplitude of each wave Pj is the full coefficient multi- 
plying the exponent e~xp(iki Z z) in the out-solution. With 
the help of the boundary conditions (8), the amplitudes 
of all the waves are expressed through the amplitude of 
the incident wave. Then after some transformations, the 
amplitude reflection coefficient, defined as the ratio of 
pressure amplitudes in the reflected and incident waves, 
can be expressed in the form 



fz-Z- iS 
fz + Z-iS' 



(27) 



Here the following notations are used. Z is a real gener- 
alization of impedance 



Z = Z g cos 9 so i, Z g = Zqx, 



(28) 



where Zq = (pos)/(p so is so i) is the ordinary impedance 
of the interface at zero frequency; \ — u/sk g is the di- 
mensionless frequency as introduced in Eq. (18); S is a 
dimcnsionlcss constant 



f (A) = k -±^ s R, 



(29) 



which is real due to Eq. (14); 

fz = .f3z/(k g f 2 z), where 
fnz =ki z (k 2z -k 3z ) + k% z (k 3z 
for n = 2, 3. 



~kiz) + k 3z (ki z - 



-hz) 
(30) 



As s so i > the transmission angles for all the waves Qi 
are less then 9 so i, so ki Z £ R and f z is a dimensionless 
real function of x an d k T /k g . 

The full transmission coefficient is = l+r_y ; the par- 
tial transmission coefficients are defined as the ratios 
of pressure amplitudes of each of the three waves in he- 
lium Pi to the pressure amplitude of the incident wave. 
They are obtained in the same way as r^: 



t; 



t. 



(ki z kj^ikiz kkz) 



(31) 



where ipi = (k iz - k +z ){k iz + k_ z ) 



and the subscripts take values {i, j, k} — {1, 2, 3}+perm. 
(perm, is for permutations). Henceforth the subscripts 
{i,j, k} in the expressions of the kind of Eq. (31) take 
the same set of values, unless stated otherwise. 

All the amplitude coefficients are complex-valued func- 
tions of frequency and incidence angle. Therefore there 
are always nontrivial phase shifts between the incident, 
reflected, and transmitted waves. 

The energy reflection and transmission coefficients are 
the normal components of energy density flux, expressed 
as fractions of the incident energy flux, that are reflected 
or transmitted into helium. The energy density flux, in 
a wave packet in the quantum fluid, as shown in [29], 
equals the average energy density multiplied by the group 
velocity. It was shown in [30], that the average energy 
density in a wave packet or plane wave in the quantum 
fluid with velocity amplitude Vi is given by the same re- 
lation as in the ordinary liquid, po|^| 2 i an d from Eq. 
(17) Vi = Pi/{poSi). The group velocity of wave i can be 
obtained from Eq. (1): 



\(kf kj)(k i 



kl) 



(32) 



Taking all of this into account and with the help of Eqs. 
(31), after some transformations we obtain the fractions 
of the normal component of the incident wave packet's en- 
ergy flux that are carried by waves of each type i = 1, 2, 3 
in helium. Those are the partial energy transmission co- 
efficients 



4-Z ki z (ki z + kj z )(kiz + fcfcz) 

(Z + f z ) 2 + S 2 kg (k iz - kj Z )(hz ~ kkz)' 



(33) 



We note that Dt* > for all i = 1, 2, 3. The full energy 
transmission coefficient can be expressed in the form 



4Zf z 



{z + f z y + 5 2 



(34) 



The energy reflection coefficient is = 



2 . Then 

from Eqs. (27) and (34) after some algebraic transforma- 
tions we can show explicitly that 

R^+D^ = 1, (35) 

and so energy is conserved when waves go through the 



8 



0.020 



0.015 



0.010 



0.005 



0.000 




10 



11 12 

#u)/k B (K) 



13 



FIG. 5. Energy dependence of the transmission coeffi- 
cients at 6 = 0, with energy in temperature units, for the 
parameterised dispersion curve (see Fig. I) and Zq = 0.01. 
The R~ roton creation probability, D2, is small. 



interface. This also verifies that the coefficients and 
R-> are the probabilities of the creation of the corre- 
sponding quasiparticles at the interface. 

Superfluid helium has a very small density and sound 
velocity, such that at the interface with a solid, the strong 
inequalities s«s OT i and po^Lpsoi hold. Then, taking into 
account Eq. (24), we have a set of small parameters 



Z < 1, Si/s sol < 1. 



(36) 



It can be shown that, due to the first condition of Eq. 
(36), in the sums of Eqs. (27), (34) and (33) the quan- 
tity Z can be neglected in comparison with the other 
terms. The second condition of Eq. (36) implies that, 
due to Eqs. (23), all the transmission angles into the 
helium are very small, as is indeed well-known for the 
interfaces between Hell and solids. Then k 2 . -C kf,k± 
and due to Eqs. (7) and (14) we have k± z ~k±. In this 
approximation we obtain f z w / = f z (x, ©so/ = 0) and 
5 ~ 5 = (k + — k-)/ik g . Then the dependance of D~* 
on the incidence angle, from (33), is factorized out and is 
reduced to the multiplier ~ cos 9 so i : 



two corresponding multipliers is reversed. So, for the 
ratio / D^ 3 , the effect is squared and we obtain 



D2 « 2V 



(39) 



Near the roton minimum, when x^Xrot, the Rr and R + 
roton branches merge, their group velocities tend to zero, 
and so do their creation probabilities D2.3 (because 
they are proportional to the energy density fluxes, which 
are proportional to the group velocities). In Eq. (38) the 
multiplier (k^—k 2 ) comes from the group velocity. In the 
same way D^ 2 ~~ ^0 near the maxon maximum x^Xmax, 
where the phonon and R~ roton branches merge. Thus 
D~2 becomes zero at both ends of the frequency interval 
in which is defined. Both the strong inequality (39) 
and asymptotic behavior of D 2 at x — > XrotiXmax (see 
Fig. 5) are the consequences of the simple relations < 
ki z < (— k 2z ) < k% z from Eqs. (20), which reflect the 
qualitative behavior of the dispersion curve for superfluid 
helium, as shown on Fig.l. 

The creation probability of R~ rotons at the interface 
is very small for all energies. It should also be noted 
that, at low temperatures, the main contribution to the 
energy flow through the interface is due to phonons of 
energies less than the roton gap (i.e. with x < Xrot), 
which are not yet taken into account. Thus twe have a 
convincing explanation why R~ rotons were not detected 
in experiments which created beams of quasiparticles in 
helium, by a solid heater, as for example [2]. 

The expressions for D~* (33), (37), and _D_>. (34) are 
written as functions of the incidence angle or k T . What 
can be measured experimentally are the energy flows as 
functions of transmission angles. If phonons be incident 
on the interface isotropically, then as the transmission 
angle for each wave is defined by Eqs. (23), the quasi- 
particles of each type are transmitted in a narrow cone 
with the cone angle twice the Of . Thus the phonons are 
injected into the helium in the widest cone and R + ro- 
tons in the narrowest cone. For the total transmission 
coefficient as function of transmission angle we obtain 



where (0)= for 6 > 6f 



(40) 



4Z g cos ( 



P + S 2 

ki z [kiz ~t~ kjz){ki z + fcfc 2 ) 

kg (kiz kj z )(ki z fcfc 2 ) 



,=0 



(37) 



The frequency dependence of the transmission factors, at 
normal incidence 6^ = 0, is shown in Fig. 5. The relative 
creation probabilities of phonons, R~, and R + rotons 
are determined by the multipliers in parenthesis in (37). 
They can be rewritten in terms of fcj, while taking care 
of the signs: k\ z and k^ z at so i = are equal to k\ and 
£13, but k 2z at 6 so i =0 is equal to (— k 2 ) (because of the 
negative group velocity of R~~ rotons) . Then for i = 2 wc 
obtain 



Do 



3C 



k 2 - ki h 



ko 



k 2 + h fc 3 + k 2 



(38) 



and k T i are the transverse components of the wave vectors 
of wave i transmitted at angle 9. 



B. Phonon or roton in the helium incident on the 
interface 



Now let us consider one of the quasiparticles of super- 
fluid helium, phonon or roton, incident on the interface. 
In terms of plane waves, a wave i of frequency lo and 
wave vector of length fej(w) is incident. The solution in 
the solid consists of only one transmitted wave, and solu- 
tion in the quantum fluid consists of one incident wave i 
and three reflected waves j = 1, 2, 3; it can be represented 
as a sum of solutions P out and from Eqs. (21) and 
(22) (see Fig. 2). The boundary conditions (8) enable us 
to express all the amplitudes through the amplitude of 
the incident wave, and thus to obtain the nine amplitude 



Both the first and second multipliers here are less than reflection coefficients r^. The coefficient nj is the ratio 
unity. In the analogous expressions for D^ 3 one of the of the pressure amplitudes of the reflected wave j to the 
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incident wave i for i, j — 1, 2, 3: 



_4h 



Ai) 
J-2z 

hz 

kiz {kf 



f® +Z-i~5 



fz + Z-i5 

- fc|) kkz/kg + Z - id 
hz f z + Z - iS 



(41) 



£ijfc.(42) 



Here the subscripts take values {i, j, k} — {1, 2, 3}+perm.; 
Eyfe is the Levi-Civita symbol, equal to 1 if {i,j, k} = 
{1,2, 3}, {2, 3,1}, or {3, 1,2} and to (-1) if {i,j,k} = 
{2, 1,3}, {1,3, 2}, or {3, 2,1}; Z = -Z g k sol z /k sol is the 
generalization of definition (28). For the incidence angles 
less than critical k\ < k 2 sol , k so[ z < and Z is given by Eq. 
(28). For greater incidence angles the new notation must 
be used, as cos# is not defined. Then fc so | 2 eC_ for the 
wave to be damped in z<0 and therefore Z = i\Z\. The 
constructions f^ nz are 

/-l = /-L/(* fl /-i). (43) 

f-L = f n z[ku,k2z,k3z]\ kit -^(-k ix ) 
for i — 2, 3; n = 2, 3. 



The amplitude coefficient of transmission for the 
incident wave of type i is 



{ki z + kj z )(kiz + kkz) 2ki z /k g 



fz + Z- i5 



(44) 



Then the energy transmission coefficient Df~ for the 
wave i can be calculated as the fraction of the energy 
of the incident wave packet that is transmitted into the 
solid. It is explicitly shown that 



Dt(x,k T ) = D^(x,k T ). 



(45) 



This important relation ensures thermodynamic equilib- 
rium between the solid and helium at equal temperatures 
on both sides of the interface. Due to Eq. (45) we can 
from now on omit the arrows in the sub- and superscripts 
of Di and D. 

The reflection coefficients for i — j are just Ru = \ru\ 
and from Eq. (41) we obtain 



(Z-iS)k g f% + f^ 



f(i) 



(Z - iS)k g f 2z + hz 



(46) 



For iy^j we have to take into account that energy flows 
for all waves are proportional to group velocities (32), 
and then from Eqs. (42) and (43) we derive 



Rij — Rji — 4kg \ ki z kj z (k^ kk)(kj fc^)] x 

2 

Z-i5 + k kz /k g 



{Z-i6)kgf 2z + hz 



(47) 



The quantity Rij is the probability of quasiparticle j be- 
ing created at the interface when quasiparticle i is in- 
cident, so the Rij can be also called "conversion coeffi- 
cients" . 

Their dependence on frequency at normal incidence is 
shown on Fig. 6. We see, in particular, that at the roton 
minimum x^Xrot, where the R~ and R + roton branches 
merge, these quasiparticles are reflected into each other 
with probability that tends to unity R23 — > 1. The same 
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effect is present for phonons and R rotons at the maxon 
maximum. 

The angular dependence of R^ is most easily analyzed 
in terms of k T instead of the three angles of incidence. 
The values of k T equal to k so i((j) or ki(uj) correspond 
to different critical angles of incidence. So, when k T G 
(0, k so i), the quantities Z and f^ nz (i.e. f nz and f^ nz 
for i = 1,2, 3) are all real, so all the waves are traveling 
waves and D 7^ 0. When k T e (k so i, fci), the wave in the 
solid is damped, Z = i\Z\ and D = 0, but f^ nz G R and all 
the waves in the helium are still reflected into each other. 
When k T € (k\ , k 2 ), the phonon wave in helium is damped 
kiz =i\kiz\ and no longer gives a traveling wave packet, 
and f^ nz also become complex. This corresponds to R ± 
rotons incident at angles greater than 21 and reflecting 
into themselves or into each other. When k T g (^2,^3) 

(i) 

the quantities f±„ z are also complex but the structure 
is different; this case corresponds to R + rotons incident 
at angles greater than 6^ and reflecting into R + rotons, 
again with probability 1. 

In all the cases, energy conservation can be explicitly 
verified but it takes different forms: 

3 

k T <k so i(u): ^2 Rij = l-Di for i — 1,2, 3; 
3 

k so i(u>) <k T <k\{uj) : Rjj = 1 fori = 1,2, 3; 
3=1 

fci(w) <k T <k 2 {uj) : R22 = R33 = 1 - ^23; 
* 2 (w)<fcr<*3(w):-Ra8 = l. (48) 

For the interface between helium and a solid, the limit 
Zq <C 1 is a good approximation, and in the Eqs. (41), 
(42) and (46), (47) Z can be neglected (in this limit- 
ing case D — > and Qf — > 0). However, the angles are 
not small anymore, as was the case for D iy and the an- 
gular dependence of the coefficients is strong. This can 
be clearly seen in Fig. 7, where the graphs of Rij and R 2 j 
are shown for hu/kB = 10K. The coefficient R 2 \ becomes 
zero at the critical angle #21 ■ The peak of R 22 and min- 
imum of i?23 correspond to angles above critical, where 
k\ z is imaginary and the damping depth of the phonon 
wave becomes roughly half of the damping depth of the 
nonlocality kernel h(r); then the imaginary part of the 
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FIG. 7. [(a) and (b)] Reflection coefficients Rij and i?2j 
for j = 1, 2, 3 as functions of incidence angle 6 at \ — 0.2871 
(hw/k B = WK). 



numerator of Eq. (47) turns to zero, and as Z is small, 
i?23 has a deep minimum. 

A more extensive analysis of the functions Rij(x, 9%) for 
the case Zq <C 1 allows us to state the following: the main 
processes near the roton minimum >Xrot+0> for all an- 
gles, are the conversion of R~ and R + rotons into each 
other and reflection of phonons into themselves; near the 
maxon maximum phonons and R~ rotons are converted 
into each other and reflection of R + rotons into them- 
selves. For phonons and R~ rotons, when the incidence 
angle becomes close to tr/2, the probabilities of reflection 
into themselves i?n j 22 tend to unity; for R + rotons this 
happens at 8 3 — > 9" 2 — 0, and at greater angles R33 — 1 
exactly. The conversion coefficients R\j j\ for j = 2,3 
are monotonically decreasing functions of the angles o f 
incidence ; Rxj becomes zero at 6j—^Tt/2 as y/n/2 — 81, 
Rjx at 6j -> 9% as ^9^-6 ]. A little above 9f x the co- 
efficients i?22 and i?33 have high sharp peaks, and i? 2 3 a 
corresponding minimum, as described above. 

Then for the case depicted in Fig. 4, the most power- 
ful beam will be always beam ii (basically because of 
phonons with energies less than Xrot)- We have shown 
that R~ rotons are hardly created by a solid heater (39), 
and the probability of R + rotons creation is also quite 
small at frequencies near Xrot (see Fig. 5) if the inci- 
dent beam consists mainly of low energy phonons. It 
was shown in [7, 8] that in a phonon beam, low energy 
phonons (1-phonons) are converted into phonons with en- 
ergy about 1QK (h-phonons). The fraction of the energy 
in the initial beam, which is converted to the h-phonons, 
can be up to 50% [8] . The conversion coefficient of these 
phonons to R~ rotons is given by R12 at fiw/fce « 10K, it 
is much greater than at the roton minimum and almost 



reaches 1/2 at normal incidence, which is more than Rxx, 
see Fig.7. 

We suggest the experimental setup depicted in Fig. 8. 
The heater injects a phonon beam, in which h-phonons 
are created. The h-phonons, incident on the solid-helium 
interface, are reflected into three beams of phonons, R~ 
rotons and R + rotons of comparable intensities (the R~ 
rotons are reflected backwards). These beams propagate 
towards the free surface of helium, and quantum evapo- 
rate atoms from it (the R~ rotons evaporate atoms back- 
ward [3]), which are then detected. Thus the energy is 
transported from the heater to the interface by phonons, 
and then to the detector by R~~ rotons along a Z-shapcd 
trajectory, with retro-reflection at the point of creation 
of R~ rotons and retro-refraction on the surface. The 
angles and fractions of the initial beam's energy, which is 
transferred to different reflected beams, are shown for the 
h-phonon part of the incident beam. The 1-phonons for 
the most part are directly reflected and are not shown. 

If the source of quasiparticles has more rotons in the 
incident beam, as the one used in [3], then beams 32 and 
23 may become also be detectable. 

The main contribution to the energy flow through 
the interface at low temperatures can be expected to 
be made by phonons below the roton gap, i.e. with 
XG (0, Xrot)- The problem of transmission through inter- 
faces by phonons with anomalous dispersion was solved 
in Refs. [30] and [20]. In the current work, the disper- 
sion relation (1) that is used is more general than the 
one used in the previous works, it is non-monotonic and 
normal below the roton gap. 

When x < Xrot, the roots k 2z ,3z are defined so that 
k\ z > and k^ z = —k 2z € C+, as shown in section 2.3. 
With these ki Z the out- and in-solutions are constructed 
(21), (22). So the amplitude coefficients are still defined 

by Eqs. (27) and (31), but the quantities /±„ z are now 
complex. The only valid reflection coefficient Ru is de- 
fined by (46) with the complex ki Z as introduced above. 
The transmission coefficient is D p ] t = 1 — R\\. It can be 
shown that in the limit of small frequencies >0, when 
the dispersion is almost linear, the expression for D p h 
approaches the standard one for linear dispersion. At 
X— > Xrot— it decreases rapidly to less than half because 
of the increasing influence of the roton waves. The curve 
-D(x) is continuous at Xrot but has a kink. 
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FIG. 8. The predicted creation of R~ rotons by h- 
phonons incident on the interface with a solid. There 
should be backward reflection and quantum evaporation 
with backward refraction. 
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IV. ENERGY FLOWS THROUGH THE 
INTERFACE 



When a phonon in the solid, of frequency u and wave 
vector k so i, is incident on the interface at angle 6 so i, 
the average energy transferred into helium is hiuD(u}, k T ), 
where D is given by (34) and k T = k sot cos 8 so i . Let the 
phonons in the solid be in thermodynamic equilibrium 
at temperature T. Then the normal component of the 
density of energy flow through the interface is (see for 
example [11]) 



Q(T) = J 



d 3 k so i 

ln !°, n T (uj)s so i cos 6 so i D, (49) 



where nx is the Bose-Einstein distribution function and 
the integration domain is the half-space k so i z > 0. The 
parts of this energy flow, that are transferred into helium 
by either phonons, or R~ , or R + rotons of helium that 
are created at the interface by the incident phonons, are 
obtained in the same way. But instead of the full coeffi- 
cient D, we now use the partial transmission coefficients 
Di . These are the corresponding creation probabilities of 
the quasiparticles. After changing the integration vari- 
ables to the arguments of Di{oj,k T ), the partial energy 
flows can be expressed in the form 



>(T) 



did 

8^ 



hui tit(oj) 



dk 2 T D % {u},k T ). (50) 



Here the upper limit by k\ corresponds to the maximum 
transmission angle of quasiparticles of type i, equal to 
Qf from (25). The quantities Qj^ for i = l,2, 3, are the 
individual contributions of phonons, R~ , and R + rotons, 
to the energy flux from the solid into helium. 

Their contributions to the energy flux, in the opposite 
direction Q*~(T), are the normal components of the en- 
ergy fluxes from helium into the solid. These are realized 
by helium quasiparticles of type i incident on the inter- 
face. The average energy transferred into the solid per in- 
cident quasiparticle of type i is, due to (45), hioDiiuj, k T ). 
Then the energy flux is derived in the same way as Eq. 
(49), with the difference that instead of s so i in the inte- 
gral we have |ttj|, because the number of quasiparticles 
incident on the interface per unit of time is proportional 
to their group velocity: 



Qt(T) = 



d 3 k 
W) 



huj nx(uj)\ui \ cos &i Di. 



(51) 



When changing the integration variables to (lj, k T ), we 
use the Jacobian determinant and we obtain explicitly 
that Qt(T) = Qr(T). 

Fig. 9a shows the ratio between the contributions to 
the energy flux through the interface of all quasiparti- 
cles above the roton gap (i.e. with huj/kg > A) Q > = 
Q1+Q2+Q3) and the contribution of phonons below the 
roton gap Q < , which is obtained using (50) with D p h- 
We see that at temperatures T < 1 K the phonons are 
dominant. However, at Tm2.5K the two contributions 
are equal, and at higher temperatures the phonons and 
rotons above the roton gap play the main role in heat 
exchange with the solid, see Fig. 9b. The contribution of 
the R + rotons to Q> increases with temperature and at 
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FIG. 9. (a) The ratio of the contributions to the en- 
ergy flow through the interface by quasiparticles above 
roton minimum, Q>, to the phonons below roton mini- 
mum, Q<, as a function of temperature (on a logarithmic 
scale), (b) The contributions of phonons, R~ rotons and 
R + rotons to the energy flow created by quasiparticles 
above the roton minimum, as functions of T 



T«3 K their contribution surpasses that of the phonons, 
see Fig. 9b. The contribution of the R~ rotons, is ap- 
proximately constant and is no greater than 6%. This is 
due to the low creation probability D 2 for all frequencies 
(39). At T = 3 K the contribution of the R~ rotons, to 
the full energy flow, is 3%. 

When there is energy flow through the interface, it 
induces the Kapitza temperature jump at the interface 
(see for example [11]). The contributions of quasiparticles 
of each type to this jump are obtained by differentiating 
Eq. (50) with respect to T. 



V. CONCLUSION 

In this work we have solved the problem of the inter- 
action of Hell quasiparticles, i.e. phonons, R~ rotons, 
and R + rotons, with the interface between helium and 
a solid. These excitations have the non-monotonic dis- 
persion curve shown in (Fig.l). The consistent solution 
of the problem has been introduced, which allows us to 
rigorously describe the simultaneous creation of the three 
types of Hell quasiparticles by any one of them, or by a 
phonon in the solid which is incident on the interface. 

When a phonon in the solid is incident on the interface, 
it is reflected with some probability and a phonon, R~ 
roton, or R + roton are created with the corresponding 
probabilities, in the helium. It is shown that the created 
Rr roton, due to its negative group velocity, is refracted 



12 



backward (Fig. 3). When some quasiparticle of helium is 
incident, all the quasiparticles, with the same frequency 
and transverse wave vectors, are created. The set of six 
critical angles as functions of frequency are introduced 
(25), (26). These separate the intervals of angles of in- 
cidence for the different quasiparticles, from which other 
quasiparticles can be created. It is shown when a phonon 
or R + roton is incident, the R~ roton is retro-reflected 
(i.e. reflected backwards), and likewise when an R~ roton 
is incident, the phonon and R + rotons are retro-reflected 
(Fig. 4). This effect is the Andreev reflection of phonons 
and rotons. 

The probabilities of creation of all quasiparticles at the 
interface, when any quasiparticle is incident, are derived 
as functions of frequency and incidence angles (33), (46), 
(47) and Figs. 5, 6, 7. It is shown that the creation proba- 
bility of an R~ roton by a phonon in the solid, and vice 
versa, is very small for all angles and frequencies (38), 
(39). This means that R~ rotons are as badly created 



by a solid heater as they are poorly detected by a solid 
bolometer. This explains the failure to detect R~ rotons 
in direct experiments until 1999 [3]. New predictions are 
made for experiments with beams of phonons and rotons 
interacting with the solid interface, and in particular, cre- 
ating R~~ rotons at the interface by a beam of high-energy 
phonons (h-phonons). 

The full energy flow through the interface, is also cal- 
culated as a function of temperature of the solid, as well 
as the individual contributions of the phonons, R + , and 
R~ rotons (50) to it, see Fig. 9. The contribution of the 
R~ rotons is shown to be very small. 
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